Unitary partitioning in general constraint preserving DAE integrators
نویسندگان
چکیده
A number of numerical algorithms have been developed for various special classes of DAEs. This paper describes a new variable step size, constraint preserving integrator for general nonlinear fully implicit higher index DAEs. Numerical implementation issues are discussed. Numerical examples illustrate the effectiveness of the new method.
منابع مشابه
Constraint Preserving Integrators for General Nonlinear Higher Index Daes
In the last few years there has been considerable research on numerical methods for diierential algebraic equations (DAEs) f(x 0 ; x; t) = 0 where f x 0 is identically singular. The index provides one measure of the singularity of a DAE. Most of the numerical analysis literature on DAEs to date has dealt with DAEs with indices no larger than three. Even in this case, the systems were often assu...
متن کاملStable Unitary Integrators for the Numerical Implementation of Continuous Unitary Transformations
The technique of continuous unitary transformations has recently been used to provide physical insight into a diverse array of quantum mechanical systems. However, the question of how to best numerically implement the flow equations has received little attention. The most immediately apparent approach, using standard Runge–Kutta numerical integration algorithms, suffers from both severe ineffic...
متن کاملSymplectic-energy-momentum preserving variational integrators
The purpose of this paper is to develop variational integrators for conservative mechanical systems that are symplectic and energy and momentum conserving. To do this, a space–time view of variational integrators is employed and time step adaptation is used to impose the constraint of conservation of energy. Criteria for the solvability of the time steps and some numerical examples are given. ©...
متن کاملA characterization of orthogonality preserving operators
In this paper, we characterize the class of orthogonality preserving operators on an infinite-dimensional Hilbert space $H$ as scalar multiples of unitary operators between $H$ and some closed subspaces of $H$. We show that any circle (centered at the origin) is the spectrum of an orthogonality preserving operator. Also, we prove that every compact normal operator is a strongly orthogo...
متن کاملFree and constrained symplectic integrators for numerical general relativity
We consider symplectic time integrators in numerical General Relativity and discuss both free and constrained evolution schemes. For free evolution of ADM-like equations we propose the use of the Störmer–Verlet method, a standard symplectic integrator which here is explicit in the computationally expensive curvature terms. For the constrained evolution we give a formulation of the evolution equ...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Mathematical and Computer Modelling
دوره 40 شماره
صفحات -
تاریخ انتشار 2004